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Abstract. In this paper we study the asymptotic behavior of the maximum magnitude of a 
complex random polynomial with i.i.d. uniformly distributed random roots on the unit circle. 
More specifically, let {rik}T=i be an infinite sequence of positive integers and let {zk}kLi be 
a sequence of i.i.d. uniform distributed random variables on the unit circle. The above pair 
of sequences determine a sequence of random polynomials Pn(z) = Yik=i ( z ~ z k)" k with 
random roots on the unit circle and their corresponding multiplicities. In this work, we show 
that subject to a certain regularity condition on the sequence {nt}'kL 1 , the log maximum 
magnitude of these polynomials scales as snI* where s% = Y^k=i n k anc ^ ^* * s a strictly 
positive random variable. 



Random polynomials are ubiquitous in several areas of mathematics and have found several 
applications in diverse fields such as random matrix theory, representation theory and chaotic 
systems (see [H [U [TJl El ED 12] ) • The geometric structure of random polynomials is of significant 
interest as well. Constructing a random polynomial from its roots is a natural construction 
which can be expected to occur in a wide range of settings. For instance, consideration of the 
asymptotic behavior of the maximum magnitude allowed the authors to obtain a lower bound 
on the minimum singular value for random Vandermonde matrices (see [13] for more details). 

In this work we continue the investigation of such polynomials, where we allow for non-constant 
multiplicity of the roots. We show that providing this sequence satisfies a simple sufficient 
condition, the limit distribution of the maximum magnitude on the unit circle is determined 
by a certain Gaussian process obtained from the Brownian bridge. Surprisingly this limit, up 
to renormalization, does not depend on the sequence itself. 

Our construction for these random polynomials is as follows. Let {^fc}^i be an infinite se- 
quence of positive integers and let {z^lfcLi be sequence of i.i.d. uniform distributed unit mag- 
nitude complex numbers. The above pair of sequences, then determine a sequence of random 
polynomials P/v(z) = Y\.k=i ( z ~ z k) Uk with roots on the unit circle and their corresponding 
exponents. Our main result relies on a construction based on the Brownian bridge and enables 
us to conclude that 
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converges weakly to a positive random variable I* where = X^fc=i n \- 
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2. Random Polynomials 

2.1. Pointwise Convergence and Lindberg Condition. Let P/v be as before and let L^(ip) 
be defined as 

N 

L N {ip) := log \P N (e^)\ 2 = J>* lo § ( 2 C 1 " cos (V> " Ok)) ) 

fc=l 

for ^ G [0, 27r]. Also let := XlfcLi n fc an d T^{ip) := Lm{iP)/sn- 

As explained in the introduction we are interested in the behaviour of the maximum magnitude 
squared of P/v(z) as N increases. Since |P/v(z)| 2 is a continuous function on the unit circle, it 
follows that there exists (p* that attains its maximum. Let Tt, be this value. For later use, let 
<3? := {ip r : r > 0} be the set of 2n times the dyadic rationals on the interval [0, 1]. Then it is 
clear that 

lim sup T N ((p r ) = Tjy. 

The case = 1 appears in connection with the asymptotic behaviour of the minimum eigen- 
value of random Vandermonde matrices (see \14\ [TH| for more details). In this special case, 
weak convergence to the normal distribution holds 

T N (iP)^ N(0,a 2 ) 

for every fixed ip where 

1 f 27T 

a 2 := — log 2 (2(1 - cosip))difj ~ 3.292. 

This is in fact a consequence of the central limit theorem since log (2(1 — cos^)) is square 
summable and Jq™ log (2(1 — cos ip)) dip = 0. 

In what follows we derive a simple sufficient condition for the asymptotic normality of the 
random variable T\r. We use the notation E (A; A) := E (XI a) where A is a random variable 
and A is a Borel set. Let A& be a sequence of independent zero mean and variance o\ random 
variables. We say that this sequence satisfies the Lindberg condition [8] if and only if 

1 N 



lim -g- V E(A 2 ; \X k \ > es N ) =0 



'N k =l 

for every e > 0. If Xk = OfcYfc then this condition becomes 



1 N / \ 

^*rMm^U . (3 ) 



lim _ 

For our purposes we only focus on the case where the random variables are i.i.d. according 
to the distribution of Y = log (2(1 — cos(2-7rf7))) and where U is uniform random variable on 
[0, 1]. The main result of this section is the following Theorem. 

Theorem 1 (Lindberg Exponent). If 

N 

E ts N 
e n k = (4) 



k=l 
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holds for every e > then 

T N W = ^l^N(0,I a ). (5) 

SN 

Proof. By definition Y = log 2(1 — cos27r?7), where U is uniform on [0, 1], from which it follows 
that Y < log 2(1 — cos n) = log 4 and hence the moment generating function exists for all t > 0. 

Additionally, since SN,n k and e are positive it follows that, 

e(y*;Y>^) <log 2 (4)P(V>^) (6) 
and applying Markov's inequality with t > we obtain 

p(y>*)<4V^ (7) 

We thus obtain the following bound for the upper Lindberg condition, 

N / \ , 0/ a\ At N 

> s>ve\ . log 2 (4)4*^ , 



4S V ~ n k ) ~ s% ^ 



< log 2 (4)4*2_]e "fc . 

k=i 

The sum on the RHS goes to as this is the condition we assumed holds and because t and e 
are arbitrary and fixed. We now turn to the lower Lindberg condition, where Y < —^f^- Let 



rik 

£,N,k '■= £SN/n k , and 9 in the interval [— 7r/3, 7r/3]. Using the following two basic inequalities, 

cos#<l-# 2 /4 (8) 

for < \6\ < vr/2 and, 

log 2 (2 (1 - cos 6) ) < log 2 (6> 2 /2) (9) 
for < \9\ < 7r/3 we obtain that, 

E(Y 2 ;Y < -£ Ntk ) < E(log 2 (# 2 /2);Y < -^,fe) 
< E(log 2 (# 2 /2);i^ ifc ) 

where ^jv.ifc = {1^1 < v / 2e^ JV > fc / 2 }. Since J \og 2 (x)dx = x log 2 (x) - 2x log(x) + 2x, we make the 
substitution = 6/V2 and determine the above expectation to be 

E(logV);M<<5) = -^=h(S) 

7TV2 

where h(S) := 5 log 2 (5) — 2<51og(<5) + 25 and (5 := e - ^-*. Rewriting this expression we obtain 

4 ^ e -es N /(2n k ) fj^L | eSN | 2 ) . 



7T \ 4n£ rife 

Summing over fc and dividing by s 2 N we obtain that 
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Covariance Kernel K(0) 




1 2 4 5 6 

Figure 1. Graph of the covariance function K(6). 



Now applying the classical Lindberg's Theorem [1] the proof is complete. □ 

The condition is easily verified to hold when = k p for p > 0. Similarly, it is easy to show 
that this condition fails if n& = 2 k . 



2.2. Finite Dimensional Limits. The result in Section 2.1 was for the marginal distribution 
of Tjv(V') f° r a fixed value ij). However, we would like to consider the weak limit for the sequence 
7jv,r where T/v,r := T^(<j) r ). It is well known that weak convergence in the sequence space M°° 
is entailed by weak convergence of the finite dimensional distributions. For this reason, it is 
important to understand the joint distribution of s such variables. We first focus on the case 
s = 2. Define the covariance function 



K(6) := (27T)" 1 / log (2 (1 - cosV>) ) log (2(1- cos (V> + 
Jo 

for 9 6 [0, 27t]. The plot of this function is shown in Figure 2.2 Given t and define, 



V N 



tT N (ipi) + uT N (ip 2 ) 



t " 



s N 



N 



m logle 



N 



(10) 



(11) 



which is a scaled sum of zero mean i.i.d. random variables with variance, 

E (V^y) = n<( (t 2 + n 2 ) a 2 + SMuirflpi - H)). 



In what follows we denote denote V = tW + uZ where W = T/v(v?i) and Z = T/v(<^2)- We 
assume that t ^ and since otherwise there is nothing to show. It is not difficult to see 
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that 

V 2 = t 2 W 2 + 2tuWZ + u 2 Z 2 

< t 2 W 2 + \tu\{W 2 + Z 2 ) + u 2 Z 2 . 



However, the triangle inequality implies that if \V\ > es^/nk then either, \W\ > es n / \2\t\nk) 
or the corresponding inequality for Z holds (or both). However, we have already demonstrated 
that 

N 



s n fc=l 

as N — > oo. As far as the terms involving Z 2 are concerned we only need to show the above in 
the case that, 

\z\ < 4T- 

that is we wish to show, 

N 



^ n fr(2P;\W\>^-AZ\< 



but the above is smaller than 

2 N 

fc=l 1 1 K 



Using Markov's inequality as in Q we see that 



es N \ ^ A „~Vz 







W > — — j— — ) < 4Pe 2| ' |n fe , 
2\t\nk> 



Applying Q and using that 6 ~ U[0, 2ir] we see that 

ts N \ : 

2|t|rifc/ ~~ 7T 



It follows that (12) tends to as iV — > oo. 

Thus by an extension of the arguments given in the proof of Theorem [T] it can be shown 
that condition Q is sufficient for the Lindberg condition to hold in respect of the random 
variables Vn,£- Therefore, (Tjv(y?i)jTjv(^2)) =^ -W(0, S) with En = E22 = c 2 and E12 = E21 = 
K(\ipi — tp2\) as an application of the Cramer- Wold device [3]. 

Clearly the above arguments go through in the case of 3 or more variables. Therefore, the 
following result holds. 

Theorem 2. Let (ipi, . . . , <p s ) be s numbers in [0, 2ir] and let (T/v(<£>i), • • ■ , T/v(y s )) be the 
corresponding random vector. Then, 

(T N (<p 1 ),...,T N ((p a ))=>N(p,'E 8 ) 

i.e. asymptotically joint normal with covariance determined by £^ = K(\(pk — ipi\) where 

K{9) := (27T)" 1 / log(2(l-cos^))log(2(l-cos(^ + 0)))#. 
J 
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Figure 2. I v for a realization of the Brownian bridge. 



3. Asymptotic Distribution of T 



N 



The following is an alternative way to construct the limit distribution of the random sequence 
T)v,r- Given a realization of the Brownian bridge W° on [0, 2tt] (which satisfies iy o (0) = 
W°(2tt) = 0). A ip shift of the Brownian bridge is defined as 



W°J6) ■, 



w°(ip + e) -w°(<p) e e [o,2-it - p], 

W°(p + 9 - 2vr) - W°(ip) 9 G [2ir — ip, 2tt]. 



In addition, define the function / : [0, 2ir] 



by 



I 



tp ■- 



2- 



W°J9) 



sml 



1 — cos ( 



-de 



for ip £ [0, 2ir]. Figure 2 shows us a realization of J^,. The following Lemma shows that / is a 
well defined quantity almost surely. 

Lemma 1. Given a realization of the Brownian bridge W° , then a.s. the following integral 
exists for all ip G [0, 2ir) 

r-2-K 



r 1 — cos ip 



< oo. 



In addition, the function is continuous. 



Proof. When <p = we write the above integral as I. The Levy global modulus of continuity 
tells us that for standard Brownian motion B on [0, 2ir) 

r v \B(t + 5)-B(t)\ 
lim limsup tt~t = 1 

< 5 ^0 0<<<27r-<5 W{d) 
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where w(8) = y 28 log j (see [TO] for a proof of this result). Since W° is by definition, 

W°{tP) = B{tP) - ^-B(2tt) 

our argument is the same no matter which value of <p is chosen because the Levy modulus 
applies to the entire sample path. We therefore set p> = 0. By definition of the Levy modulus, 
there exists 82 > almost surely such that 

\B(t + 8) — B(t)\ 



w(8) 



< 2 



for all < 8 < 82- Therefore, 



a(8) := \W°(tp + 8)- W°(tf;)\ < 2w{8) + 



2tt 



(13) 



We may therefore split the integral as, 



2n-8 2 



82 



W°(ip) 



simp 
1 — cos Ip 



82 



dtp + / W°{tp) 



simp 
1 — cos tp 



dtp + 



2tt 



2tt-5 2 



\v°ty) smip M . (14) 

1 — cos ip 



The first integral is finite being the integral of a continuous function over the interval [82, 2tt— 82]- 
We may further suppose that 82 has been chosen so that \ip | < 4 for < tp < 82 with 

the corresponding inequality in a similar neighbourhood of 2ir. By choice of 82 we obtain that 



V(*) T ** 7 # 

1 — cos tp 



< 



6 2 



= 0(8\") 



for sufficiently small 82- The same argument applies to the last integral. Since 10(^2) gives 
a uniform bound the result holds for all 92 6 [0, 2tt). Continuity in ip follows by a similar 
argument, 



< 



2n-8 



cos ip 



(15) 



W£(tp)\ 







sin^ 



1 



+ 
+ 



2n \w^)\ 

2n-5 



cos ip 
simp 



dtp 







1 — cos ip 
sin ^ 



2 77 



2tt-8 



1 



cos ip 
simp 



1 — cos ip 



dtp. 



Provided that < 8 < 82, the tail integrals are all at most 0{8 l /' i ) as before. We bound the 
first integral by two positive integrals, to obtain 



27T-5 



Wr) ~ 



svntp 
— cos ip 



dtp 



< 2sup|W°(^)-W|(V)| 



smtp 



dtp 



l s (1- cost/0 

< 2 sup I W^tP) - W^tP) I [ log (1 - cos rP) ] ] 

< 6a(8) (log 2 - log(l - cos 8)) 



s 
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provided \tp — (p\ < 5. Finally, since 

a(S) (log 2 - log (l - cos 5) ) -»• 
as 5 — > we finish the proof. □ 



Let <3? = {(p r : r > 0} be the sequence described in Section 2.1 and let I = {/ r }^ be the 
sequence defined as I r := I Vr . Since the function I v is continuous on the interval [0, 27r] there 
exists a value ip* which determines the maximum value of 1^, which we denote by I*. Since $ 
is dense on the unit circle it follows that 

I* := sup {I r : r £ N} (16) 

and its distribution is determined via the infinite sequence I r . Let the sequence of random 
variables TV = {Tjv>}^o- We now derive one more Lemma for use later on. 

Lemma 2. Let Y be a function in D[Q, 2tt]. Then Y is Lebesgue measurable, and its integral 
exists, 

/ Y(s)ds < oo. (17) 
Jo 

Furthermore, let Y n be a sequence of functions in D[0, 2n] such that Y n —> Y in D (i.e. with 
respect to the Skorohod topology) then 

2n r-2ir 

Y n (s)ds^ / Y(s)ds. (18) 



Proof. The existence of the integral follows from Lemma 1, page 110 of [3 J and the subsequent 
discussion which shows that functions in D on a closed bounded interval are both Lebesgue 
measurable and bounded. The former follows from the fact that they can be uniformly approx- 
imated by simple functions, a direct consequence of Lemma 1 and the latter also. 

Convergence follows from the Lebesgue dominated convergence theorem. This holds since the 
sequence Y n is uniformly bounded, by a constant so the sequence is dominated. Second Y is 
continuous a.e. with pointwise convergence holding at points of continuity, as a consequence of 
convergence in D see [3]. □ 

We now proceed to prove the following Theorem. 

Theorem 3. The sequence T^v converges in distribution to the sequence I 

T N I (19) 

as N — ^ oo. 

Proof. In order to do so we use Theorem 4.2 of [3]. Suppose that there is a metric space S 
with metric po and sequences Tn,€, I e arid T^r all lying in S such that the following conditions 
hold, 

T N , t I e (20) 

L => I 
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together with the further condition that given arbitrary r] > 0, 



limlimsupP(/5o(TV ie ,TV) > rj 



0. 



(21) 



Then it holds that Tjy I. First, we define I e using a realization of the Brownian bridge as 
follows, 

/■27r-e 

I 2-rr-e 



I*,, ~ K,Wlog2(l-cos. 



r 1 — COS ip 



The definition of the other sequence is more involved and so we defer it for a moment. We 
have shown that the limit integrals exist a.s. and so we only need to show that the first term 
converge to 0. Since log (2 (1 — cos ip) ) = O(loge) when e is small and in a neighbourhood of 
and 2tt we may invoke the Levy modulus of continuity, wrapped around at 2tt to obtain that 
this term is 

0(logea(e)) -> 0. 
Hence, coordinate convergence of the integrals holds so that 



Ir.e + 



2ir-e 



^ 1 — cos ip 







and it follows that I e => I as e — > 0, since the sign of the integral is immaterial. We have thus 
demonstrated the second condition of (20). Next, we proceed by rewriting T^{^p r ) in terms of 
the empirical distribution function Fjsr : [0, 2ir] — > [0, 1] determined by 

#{e q -.o<e q <^} 



F N {i>) :-- 



By definition of F^(ip) and the Lebesgue-Stieljes integral we see that 



T N (<p r ) 



2tt 



sn log(2(l-cos(ip r -<ip)))dF N (ip) 



2- 



s N log (2(1 - cosip))dF Njipr (ip) 



where the change of variables, ip = ip — if has been made. For ip G [0, 2tt) we define F^ )lf {ip) as 
the "cycled" empirical distribution function of F^ by 



#{<p<0 q «p+if>} 

SN 



F N , v (2ir) + 



#{0<6> 9 <V-2 7 r+ ¥ >} 



if ip < ip < 2ir — ip, 
if 2tt - ip < ip < 27T. 



2tt 



To define the sequence T^ )€ (<p r ) we split the integral into two parts as in f e n e and J e + f 2 * 
and then use integration by parts on the first part, which yields the expression, 



TN,e(¥r 



'SN 



Fn«M) 



2vr 



~1 2tt-£> 



log 2(1 — cos ip) 



(22) 



2lT-t 



simp 



2ir J (1 — cos ip) 



dip. 
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For later use we make the definition, 



This is not quite equal to the original sum, since 

r>2vr 



/ log (2 (1 — cos 
Jo 



so that the ifi terms do not give but rather cancel with fi 6 to be defined in a moment. The 
remainder we express it as a sum, noting that we must include the mean, which is by symmetry, 

Ve -=~ [ log (2(1 -cos 2 (eloge-e + o(e)). (23) 
Zir Jq 7t 

Define S e (tp) := {9 q : 9 q G [(p — e, ip + e]} and hence the sum can be written as 

Z N ^ip r ):=-= V log (2(1 - cos(^ r - ^))) - -=He- (24) 

\/SN \/SN 

Denote the corresponding sequence as Zat j(E . Taking expectations we thus find that 

E(Z N 4ip r )) = -L / log(2(l-cosV))#-^M e = 

is a sequence of random variables with mean. We finally write, 

Tat = Tjy.e + Zjv, e . (25) 



We now proceed to demonstrate the first condition of (20), namely that, T^,e =^ le- The 
random variable T/v, e (<£v) is a functional of an empirical distribution and therefore of a process 
lying in D[0, 2tt]. Define the random sequence J € defined for / £ D[0, 2tt] and f(0) = f(2ir) = 
with the component term, 

fZ-ir-e sin?/; / \ 27T-e 

JeAf) = j UM (1 _ cos ^ # - [UM log (2(1 - cos v))) e • (26) 

It is well known that Wn,q =>- W° in D, which implies that Wjv )¥>r W° r as N — > oo for all 
r. The result follows on showing that J e defines a measurable mapping J e : D[0,2n] — > M 00 in 
£[0,271-]. Since 

Je,r(W N ) = T N)e (if r ), 

we may therefore apply Theorem 5.1, Corollary 1 of [3] which states that if Wjv M^ then 
■^(W/v) =^ J t (W°), (and hence Tjv, e =>■ I € ) provided that we verify 

W° G Dj.) = 0. (27) 

To deal with the measurability question we first observe that the coordinate maps are mea- 
surable and since smip/(l — cos-0) is continuous in [e, 2ir — e], it follows by Lemma that 
J £iT . is measurable for each r and hence so is the sequence mapping J e . Again by Lemma]2]the 
sequence of integrals convergences with respect to pq. This leaves only the final term. However, 
since the limit W° is almost surely continuous it follows that 

UM -»• W° r {e) 
/^(2vr-e) -»■ W° r (2n-e) 
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for each r if / — > W° in D[0,27r]. Thus the corresponding sequence converges with respect to 
po also and so (27) holds. The proof of the first condition is concluded. 



It remains to demonstrate ( |2l| ). Here we use the union bound and Chebyshev's inequality. 
This is because the various ZN^((p r ) in the sequences are dependent, as they are determined 
via the same 9 q . Nevertheless they are of course themselves the sum of i.i.d. random variables. 
In determining the variance, we may work with tp r = without loss of generality. The variance 
of one of the i.i.d. summands in (24) is determined as 

o 2 r 

2W 



log 2 (2(1 - cos tp))dip -/■*?< 



oo. 



(28) 



Since for small e > we have log (2(1— cos if))) = 0(2 log ip)+o(ip) the integral is a e = 0(e log 2 e) 
as the integral of log 2 x is xlog 2 x — 2xlogx + 2x. It follows that of — > as e — > 0, which is 
the variance of the entire sum by independence and as it has been scaled. 



Now fix r] > 0. By definition of po an d from (25) we obtain that, 

Po(Tn,6, Tat) 



^2 iznAp^i 2 - r 



r=0 



+ \ZN,e{fr)\ 



Let R n be such that YlT=R n +i 2 r < f?/2. Now we apply the union bound to the remaining 



Rrj + 1 summands to obtain that 



R 



-ft) 



r=0 



+ \ZN,t{Vr)[ 



< 



< 



5>(|Z* >e (p r )|2- 

r=0 
R v 



> 



2 (^ + 1) 



(29) 



r=0 



?? 2 2 2r 



< <7. 



2 16(i^ + l)^ 



Hence, 



limSUpP(p (T iVie ,T A r) > 7]) < 



3rj 2 

16(R V + l) 2 a 2 
3rj 2 



and the RHS goes to as e to 0, for each rj > 0. Hence we obtain (21) as required. Therefore, 
we have verified all conditions and Theorem [3] is proved. □ 



It is rather easy to see that I* > almost surely. For instance, as in the proof of Lemma [T] it 
can be shown that, 

r2n i-2-k 











svntp 
1 — COS if) 



dip dip < oo. 



It then follows from Fubini's theorem that, 



2rr 



l<pd(p = 
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0.5r 
0.4; 

0.3 1 

0.2; 

0.1- 



Prob. distibution function of I* 




III. I. 



3 4 p 5 



Figure 3. Histogram of the pdf of I*. 

since W°(ip)dip = 0. But I v is almost surely continuous and hence I* = if and only if 
Iv> = for every ip. It thus follows that I* > almost surely, as required. 



Theorem 4. Given r > 0, 



liminfPIT^ > r) > P[ J* > r 



AT->-oo 

with equality if t is a continuity point for the random variable I* . Moreover, 

— log max { | P;v ( z ) 1 2 

SAT I 



1 



r. 



It therefore follows that, 

max{|P^(z)| 2 : \z\ = l} « e V"i+-+^^ 
for iV sufficiently large. A sample histogram for /* is shown in Figure [3| 



(30) 



4. Numerical Results 

In this Section, we present some simulations of our results. In Figure [4] we show the logarithm 
squared magnitude for a random polynomial with N = 500 and constant sequence n& = 1. Here 
the maximum value is ~ 10 60 and occurs near ip = 0.3. Our final plots show the logarithm 
of the maximum magnitude as a function of the degree for the sequences n& = 1 and the 
sequence = k with 100 realizations per degree. The blue curves are (n\ + . . . + n^) 1 / 2 and 
5(n 2 + . . . + n 2 ^) 1 / 2 respectively for both cases. 
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Figure 5. Logarithm of the maximum magnitude as a function of the degree 
for the sequences nt = k (top) and the sequence = 1 (bottom) and 100 
realizations per time. The blue curves are (n\ + . . . + nj^) 1 / 2 and 5(nf + . . . + 
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) 1 / 2 for both cases. 



